y<-c(2, 6, 10, 8, 12, 9, 28, 28, 35, 32, 46, 47,
          50, 61, 99, 95, 150, 143, 197, 159,
         204, 168, 196, 194, 210, 180, 277, 181,
         327, 276, 365, 300, 356, 304, 307, 386,
         331, 358, 415, 374, 412, 358, 416, 414,
         496)
x<-1:45
qrt<-gl(4,1,length=45) 
aids<-data.frame(y,x,qrt)
rm(y,x,qrt)